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ABSTRACT 

It is anticipated that the large sky areas covered by planned wide-field weak lensing surveys will reduce statistical errors to such 
an extent that systematic errors will instead become the dominant source of uncertainty. It is therefore crucial to devise numerical 
methods to measure galaxy shapes with the least possible systematic errors. We present a simple "forward deconvolution" method, 
gfit, to measure galaxy shapes given telescope and atmospheric smearings, in the presence of noise. The method consists in fitting 
a single 2D elliptical Sersic profile to the data, convolved with the point spread function. We applied gfit to the data proposed in 
the GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Galaxy Challenge. In spite of its simplicity, gfit obtained the lowest 
additive bias ( V5l = 0.057 x 10~ 4 ) on the shear power spectrum among twelve different methods and the second lowest multiplicative 
bias (A4/2 = 0.583 x 10~ 2 ). It remains that gfit is a fitting method and is therefore affected by noise bias. However, the simplicity 
of the underlying galaxy model combined with the use of an efficient customized minimization algorithm allow very competitive 
performances, at least on the GREAT 10 data, for a relatively low computing time. 
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1. Introduction 



Weak gravitational lensing (e.g., Bartelmann & Schneider|20TjT 



|Hoekstra & Ja in 2008 ), whereby the gravitational bending of 
light by structures in the Universe slightly distorts images of dis- 
tant galaxies, is now recognized as a powerful means to study the 
history of the Universe and probe the mysterious nature of the 
dark matter and dark energy (Munshi et al. 2008| Huterer|2010| l. 



Since the first detection of weak lensing (|Maoli et al. 



Bacon et al.|2000||Kaiser et al.|2000[|Van Waerbeke et al. 



2001 



2000 



Witt man et al.|[2000[ ), a number of methods have been devised 
and implemented to tackle the inverse problem of recovering 
the lensing signature from observed, distorted galaxy images 
(|Kaiser et al.|| 1993| Luppino & Kaiser||1997| |Hoekstra et al.| 



1998 , Bernstein & Jarvis 2002; Hirata & Seljak 2003; Refregier 
& Bacon|2003||Heymans et al.|2006||Massey et al.|2007[|MiIler1 
et al.||2007[ |Kitching et al.||2008[ |Bridle et al.||2010| |J^tching| 
et al.|2012) . 

We describe in this paper gfit, a simple shear measure- 
ment method that nevertheless obtained good results in the lat- 
est GRavitational l Ensing Accuracy Testing 201 (GREAT 10) 
Galaxy challenge (Kit ching et al.|[2011| |2012| l. Galaxies are 
assumed to be well modeled by a seven-parameter, single- 
component elliptical Sersic profile. The shape measurement al- 
gorithm essentially consist in iteratively shearing and convolving 
the galaxy model until a sufficiently close match with the ob- 
served galaxy is reached. Instead of an out-the-box minimizer, 
we employ a custom-developed minimizer well suited to fitting 
faint and noisy images like those frequently found in weak lens- 
ing. 

The paper is structured as follows. We provide in Sect. [2] a 
description of the underlying principles, galaxy model and shape 
measurement algorithm of gfit. We continue in Sect.[3]with a pre- 



sentation of the pipeline we used to participate in the GREAT10 
Galaxy challenge and follow with a analysis of the gfit results in 
Sect.g] We conclude in Sect. [5] 

2. The gfit shear measurement method 

2.1. The shear measurement problem 

According to the theory of weak gravitational lensing, the light 
emitted by a galaxy is slightly deviated by the foreground grav- 
itational field, an effect that can be modeled to first order as the 
combination of two effects, the convergence k and the shear y, 
that describe how light bundles emitted by a source are distorted 
by a potential well. The convergence models the magnification 
effect whereby the galaxy image see its apparent size increased 
without altering its shape, whereas the shear describes a stretch- 
ing effect where only ellipticity is altered, not size. 

All so-called "shear measurement methods" attempt to re- 
construct the reduced shear g = y/(l - k — y) which is approx- 
imately equals to the shear y in the weak gravitational limit, 
where k <sc 1 and y« 1. 

The lensing effect is very subtle, however, and requires mea- 
suring the shapes of thousands of faint galaxies. Moreover, be- 
fore they reach the observer, the apparent galaxy images undergo 
a number of additional distortions, unrelated to lensing, that fur- 
ther complicates that task, mainly: 

- The convolution of the images by the instrumental and/or at- 
mospheric point spread function (PSF) that flattens and cir- 
cularizes the galaxy light profile. 

- The Gaussian and Poisson noise introduced by the surround- 
ing sky emissions and the detecting device. 

- The pixelation effect caused by the integration of light falling 
on the detector pixels. 
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The traditional approach for estimating the shear is to mea- 
sure the deviation from circularity of a large number of galaxy 
shapes. But this technique assumes the shear remains con- 
stant across the field of view, which is generally not the case. 
Moreover, accurate shear measurement must also account for a 
spatially varying PSF that must be interpolated at the positions 
of the galaxies. 

A shear measurement pipeline must overcome all the above 
difficulties, typically going through the following steps: 

1. PSF correction, whose goal is to restore the shape a galaxy 
had before being convolved with the PSF. If the spatial 
variation of the PSF over the field of view is significant, this 
step also requires interpolating the PSF to the position of the 
galaxies in the sky. 

2. Shape measurement, that is, the estimation of the galaxy 
shape after it has been altered by the cosmic shear but before 
PSF convolution and other subsequent distortions. In this 
paper, we call the corresponding shape the sheared galaxy 
shape. The real galaxy ellipticity prior to gravitational lens- 
ing (i.e. unsheared) is referred to as the intrinsic ellipticity. 

3. Shear measurement, that is, the task of extracting the shear 
signal from the sheared galaxy shapes estimated in the pre- 
vious step. A spatially-varying shear field is commonly de- 
scribed as a power spectrum or a correlation function. 

Additional steps may also be performed to correct the images 
from the effects of pixelation and noise. 



2.2. Shear measurement with gfit 

The gfit "galaxy fitting" method grew from a prototype initially 
developed by Stephane Paulin-Henrikson on the occasion of the 
GREAT08 challenge ( |Bridle et al.|2008| ), where it obtained the 
third and fifth ranks on images with high and low signal to noise 
ratios respectively (Brid le et al.|2010 1. 

The GREAT10 gfit code was subsequently made more 



generic in order to satisfy the more demanding requirements of 
the Galaxy challenge. It was also enhanced in several aspects 
that we describe in subsequent sections of this paper. 

It is also worth mentioning that a wavelet-based denoising 
algorithm described in Nurbaeva et al. ( 201 l| l was also exper- 
imented in the Galaxy challenge and proved quite successful. 
This denoising scheme is presented in Sect. |2.7| 



The overall shape measurement procedure is the following: 

1 . Application of the denoising algorithm on the galaxy and/or 
PSF images [optional] 

2. Estimation of the galaxy and PSF centroids in all images 

3. Application of the PSF correction and shape measurement 
algorithm 

4. Generation of the ellipticity catalogs 

5. Production of various statistics and plots for analysis [op- 
tional] 

We detail in the next sections the shape measurement algo- 
rithm along with its underlying models and components. 



2.3. Modeling the galaxies 

gfit is fundamentally a model-fitting method. We describe here 
the model used to represent galaxies and cover the fitting-related 



aspects in Sects. |2~4| and |2.5| 

Galaxies are assumed to have a surface brightness distribu- 
tion well described by an elliptical Sersic function ( Sersic 1968| l, 
defined by: 



/(£ n, re) = I sky + I exp [ - b„ (-) "] 
where: 



(i) 



/ S jh, represents the sky brightness 
Io is the central surface brightness of the galaxy 
n denotes the Sersic index that determines the degree of cur- 
vature of the profile. A small value of n leads to a less cen- 
trally concentrated profile and a shallower logarithmic slope 
at small radii. 

The scale radius parameter r e is defined as the effective ra- 
dius encircling half of the total light of the profile (e.g. Ciotti 
[T99Tj|Trujillo et al.|2004[|Graham & Driver|2005) . 
The factor b n (e.g. Ciotti & Bertin|1999 1 arises from the def- 
inition of r e and is related to the Sersic index n through the 
equation F(2n) = 2y{2n,b n ) where F and y functions are 
respectively the complete and incomplete gamma functions 
( [Abramowitz & Stegun|1965| l. 
The parameter £ in Eq. ([IJl is defined as 



Kx>-x c ) 2 



(y'-yc) 2 

q 2 



(2) 



and denotes the distance from the centroid (x c ,y c ) of the 
galaxy to a point on an elliptical isophote at spatial coor- 
dinate (x',y r ) 
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obtained after counterclockwise rotation through an angle <p 
with respect to the (0, x) axis. 
- The quantity q in expression (|2]) of £ is the ratio of the semi- 
minor axis b to the semi-major axis a of the isophote el- 
lipse. It is related to the complex ellipticity e = (61,62) 
of the galaxy through q = b/a = (1 - |e|)/(l + |e[) with 
|e| cos 2<p and e2 
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|e| sin 2(f>. 



This model was initially chosen for its simplicity and has the 
other merit of being relatively easy to fit. A galaxy represented 
as the sum of a bulge and a disc would be more realistic, but it 
is not clear whether such a model would prove more accurate on 
weak lensing images and worth the additional complexity and 
computational cost, given degeneracies between parameters. As 
regards gfit, the GREAT10 results does not provide a definitive 
answer (see Sect. [4]). 

2.4. Galaxy shape measurement and PSF correction 

We describe in this section the shape measurement algorithm of 
gfit. It is based on iterative fitting of observed galaxies to the 
galaxy model described in Sect. 2.3 The basic assumptions of 
the algorithm are the following: 

- All galaxy and PSF fields have been reduced and the ob- 
jects they contained assumed available in the form of square 
postage stamps in FITS format. We denote by O a galaxy 
field image of the GREAT 10 Challenge and by o any of its 
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10,000 galaxy postage stamps. The GREAT 10 PSF field im- 
age corresponding to O is labeled as P and any of its individ- 
ual PSF kernel as p. 

- The O and P images have been denoised beforehand if re- 
quested. The denoising tool described in Sect. 2.7 is the de- 
fault choice. 

- The centroids of all galaxies and PSF to be processed 
have been estimated. A tool based on SExtractor (Bertin 



& Arnouts 1996 1 has been written for this purpose (see 



Sect.[Z6|) 

- The PSF field at the position of each galaxy is known, so that 
we can always find the PSF p that matches galaxy o. 

The fitting algorithm itself is summarized below: 

- For each galaxy object o in observed galaxy field O 

o Extract a square stamp cutout of a given dimension Nq x 

Nc around the centroid (x c ,y c ) of o. 
o Remove the sky background from o after having esti- 



mated its variance cr 2 ky . 



Estimate the galaxy noise variance a 2 of o. 

Select guess parameters for the 7 mod el pa rameters 



2.3 



{/o, (x c ,y c ), (ci,«2). n, r c } described in Sect, 
o Construct an initial galaxy postage stamp g based on the 

Sersic profile having these parameters, 
o Select the PSF p that matches the spatial coordinates of 

o in O and extract a cutout of dimension Np x Np about 

the centroid. Make sure the PSF is normalized and has 

the sky background removed, 
o Iteratively vary the model parameters with the objective 

of minimizing the residuals between g and o. 

At each step: 

Convolve the model galaxy g with the observed PSF 
p and compute a new galaxy estimate k = p ★ g 
=> Compute the residuals between o and k using the 
chi-squared statistics 

^=zz^ — 

n m 

where o„ „, and represent the pixel value at po- 
sition (n, m) in o and p respectively. Similarly, cr n m 
denote the standard deviation of the noise associated 
with each pixel at (n, m). 
=> If the minimum x 2 nas n °t been reached, select a 
new set of parameters and construct the correspond- 
ing model galaxy g. Otherwise, exit the minimiza- 
tion loop. 

o At the end of the iteration cycle, the algorithm yields: 
=> An estimate for all 7 model parameters and for the 

complex ellipticity (<?i, ej) in particular. 
=> A postage stamp for the best-fitted intrinsic galaxy 

postage stamp g. 
=> A postage stamp for the best-fitted convolved galaxy 

k. 

o The same procedure is followed for the next galaxy o in 
O, until all galaxies have been processed. 

- Once the shape of all galaxies have been measured, pro- 
duce a catalog containing the fitted parameters at the posi- 
tions of the galaxies o in O and optionally, additional statis- 
tics and plots. Quantities such as the ellipticity modulus |e|, 
the position angle <p or the minor-to-major axis ratio q can 
be respectively derived from (e\,e-i) using |e| = V e i 2 + ei 2 , 
<p = \ arctan(e 2 /ei) and q = (1 - |e|)/(l + |e|). 



The algorithm just described is not dependent on the galaxy 
or the PSF model. It would remain unchanged, for instance, if the 
centroid of the PSF itself was taken into account during fitting 
or if a bulge and a disc were incorporated in the galaxy model. 
The quality of the whole shape measurement procedure can thus 
be increased by improving the models. 

The shape measurement algorithm has its own strengths and 
weaknesses, outlined below: 

- Strengths: 

o Simplicity. 

o PSF deconvolution and galaxy model estimation are per- 
formed simultaneously. 

o The deconvolution process does not involve any ma- 
trix inversion and its side effects (numerical instabilities, 
presence of artifacts, noise amplification, etc.). 

o The intrinsic and best-fitted modeled galaxies are ob- 
tained as a by-product of the algorithm, in addition to 
the estimated model parameters. 

- Weaknesses: 

o A bias in introduced through the choice of a galaxy 
model, which is necessarily imperfect. 

o The choice of the initial guess parameters of the galaxy 
model are more or less arbitrary and may also influence 
the final model parameter estimates: 

o The algorithm relies heavily on the accuracy and robust- 
ness of the minimizer. 

o The methods is sensitive to noise bias 

o No estimate of errors is available 



2.5. The gfit minimizer 

The gfit galaxy model expressed by Eq. ([T]i varies linearly with 
the parameter /(>, but non-linearly with the remaining parameters 
{(^ c ,y c ), {e\ , ei), «, r e J. This requires the use of a non-linear min- 
imization algorithm over a seven-dimensional parameter space. 

Unlike linear minimizations schemes that only involve a ma- 
trix inversion, non-linear optimization requires iterating over the 
parameter space to find the minimum value of the objective func- 
tion, which is in our case the x 2 of me residuals between ob- 
served and estimated images. That minimum is not necessarily 
the absolute minimum of the x 1 function but the most relevant 
from the point of view of the physics of the problem. In our case, 
the minimum should coincide with the Sersic model parameters 
that best fit the galaxy shape. 

A good minimizer is essential to any galaxy model fitting 
algorithm, but finding such a minimum in a reliable manner can 
prove tricky for a number of reasons: 

- Parameter degeneracy, where different combinations of the 
parameters yield similar^ 2 values. 

- Errors related to noise and undersampling may degrade the 
accuracy of the fit in several ways: (1) the distortion of the 
galaxy image may be such that the minimizer fits a wrong 
shape, even if it does it accurately. (2) The minimizer may 
be fooled by a "false" local minimum and converge toward 
wrong fitted parameter values. (3) The minimizer may fail to 
converge altogether if it cannot reconcile the model with the 
observed image. 

- The choice of initial guess values for the parameters can in- 
fluence the outcome of the minimization algorithm: depend- 
ing on the starting location on the parameter hypersurface, 
the algorithm may tend to follow a different path and con- 
verge to a different minimum. 
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In this regards, the choice of the galaxy model used by gfit 
leads to a number of challenges: 

- The total flux Iq is degenerate with the n and r e parameters, 
so that, for instance, in Eq. ([T}, a low Iq flux may be compen- 
sated by a higher exponential function of n, re and (x c ,y c ), 
modifying the shape of the fitted galaxy profile. Similarly, an 
error in the estimation of the sky brightness I s k y may "drive" 
the minimizer toward a wrong combination of the remaining 
parameters. This is why we estimate 1^ separately in the 
algorithm described in Sect. |2.4| 

- gfit attempts to limit the effect of noise by applying the de- 



noising method described in Sect. 2.7 As regards pixeliza- 
tion, gfit can optionally construct the galaxy model in higher 
resolution and rebin the pixels before fitting. 
- To estimate the initial guess values for fitting a particular 
galaxy, gfit can be setup to either accept default values or use 
estimates from SExtractor. 

Several families of optimization algorithms were ex- 
perimented on GREAT08 and GREAT 10 images: simplex 
(Nelder-Mead downhill), gradient descent {Powell), Newton & 
quasi-Newton (Newton-CG, BFGS) and Levenberg-Marquardt 
(LVM). Descriptions of these algorithms can be found in (e.g. 
Levenberg| 1 944| |Marquardt| 1963] |Powell| 1 964| ~ 



Nelder & Mead 



1965 1 |Zhu et al.|1997[|Nocedal & Wright|1999 



2006) 



Bonnans et al. 



None of these methods proved entirely satisfactory, either 
failing to converge or yielding insufficient accuracy, especially 
on low signal to noise ratio (S/N) images. The LVM implemen- 
tation from the SciPy library ( |Jones et al.||2001^| |Corrmiuiiity 
|2010| that we used was the fastest and the most accurate. For 
these reasons we used it in the GREAT08 version of gfit, but that 
implementation of LVM: 

- Failed to converge in about 5-10% of the time on GREAT08 
"real noise blind" images. 

- Was occasionally tricked by "false" local minima, producing 
the smallest residuals but with unphysical Sersic parameters 
or ellipticities. 

- Required good estimates of guess parameters in order to con- 
verge towards the right minimum. 

We also experimented a parameter estimation scheme based on 
a Bayesian approach and implemented using the pymc Markov 
Chain Monte Carlo (MCMC) library of |Patil et alT] ( |20T0] >, but 
that method produced less accurate estimates while being less 
computationally efficient. 

In an attempt to overcome these issues, we eventually de- 
cided to implement a custom minimizer, better suited to fitting 
noisy, pixelized galaxy images than vanilla minimization algo- 
rithms. We found that a scheme based on an adaptive cyclic 
coordinate descent algorithm (CCD) was able to produce more 
accurate estimates while at the same time being more robust: 

- Able to better cope with degenerate and correlated parame- 
ters: LVM has difficulties with the {/o, n, r c J degeneracy and 
coupling and the steepest descent algorithm used in LVM 
occasionally jumps without precaution to a minimum value 
that may be the smallest but not the most appropriate one. To 
avoid this, the CCD algorithm ensures that these parameters 
are carefully varies in a "round-robin" manner at the begin- 
ning of the fitting process, where the amplitudes in variation 
are the greatest. This scheme is also much more tolerant with 
regards to initial guess values and converges reliably, making 
it more robust overall. 



- More resilient to noise and on average more accurate than 
LVM, especially on low S/N images. 

The CCD algorithm proved suitable for fitting without any single 
failure the huge number of GREAT10 galaxies. It nevertheless 
has a number of drawbacks, namely: 

- The convergence rate is lower than that of LVM, resulting 
in a greater number of function evaluation. It is thus much 
slower than LVM. 

- Its efficiency decreases rapidly with the number of parame- 
ters. The algorithm performance is also influenced by factors 
such as the specific stopping conditions chosen or the range 
if iteration step sizes specified for the parameters. 

The current gfit implementation can be configured to use ei- 
ther the CCD, LVM or MCMC-based minimizer. 

2.6. Centroid estimation 

gfit does not assume objects to be correctly centered within their 
postage stamps and accurate estimates for galaxy and PSF ob- 
jects are required for two main reasons: 

- gfit does not necessarily use the whole postage stamps for the 
galaxy and the PSF to save computational time and reduce 
noise: the corresponding postages stamps are cut out to a 
smaller dimension (e.g. 24x24 instead of e.g., 48 x48 around 
the estimated centroid. 

- The coordinates of the galaxy and PSF centroids are pro- 
vided to the minimizer as initial guess values before the fit- 
ting cycle can begin. Accurate centroids are especially im- 
portant if the LVM minimization algorithm is used (CCD is 
much more tolerant in this respect). 

gfit relies on centroid estimates obtained from the SExtractor 
tool (Bert in & ArnoutsJl 996). A catalog is generated with cen- 
troid information and additional data such as flux, ellipticities 
and position angles, that can be optionally used to set guess pa- 
rameter values. 



2.7. Denoising 

Correcting astronomical images from the effect of noise has al- 
ways been a challenging task. This is particularly true for galaxy 
images captured for weak lensing analysis because noise not 
only degrades the overall quality of these images but also al- 
ters the shapes of the galaxies. This causes serious difficulties to 
all existing shear measurement schemes that found their shear 
extraction algorithms on the accurate measurement of galaxy 
shapes. 

The challenge is then to correct galaxy images from 
noise without compromising the shear signal they encode. 
Unfortunately, popular denoising algorithms based on median 



filtering (Arce 


2005 


Arias-Castro & Donoho 2009), Wiener fil- 


tering (Wiener 


1949 


Khireddine et al. 2007 1 or discrete wavelet 



are not shape preserving and do not meet that requirement. 

By default, gfit uses DWT-Wiener, a shape-preserving de- 
noising technique combining DWT and Wiener filtering devel- 



oped at the laboratory of astrophysics of EPFL by Nurbaeva 
|et al.|p0TT) . That algorithm has been experimented during the 
GREAT10 Galaxy challenge and was able to significantly im- 
prove the quality factors of all the shear measurement methods 
from EPFL that participated in the Galaxy challenge (Kitching 
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et al. 2012|l: that w as the case for gfit but also for MegaLUT 
( [Tewes et al.|2012| > and TVNN ( |Nurbaeva et al.|in prep!) 



Interestingly, denoising improves the shape measurement of 
the three algorithms even though they are fundamentally differ- 
ent from each other. 

Beyond shape-preservation, another advantage of the 
DWT-Wiener algorithm lies in its ability to denoise "in one go" 
images containing a great number of objects, without having to 
individually process each object in turn. In GREAT10, for in- 
stance, DWT-Wiener was directly applied to images containing 
100 x 100 PSF or galaxy postage stamps. 



SET 01 








Galaxy 
Image 001 

J 



GALAXY 
CHALLENGE 



Galaxy 
Image 200 



Challenge contains 24 sets 



Each Set contains 200 galaxy images, each with 
10000 galaxy postage stamps. Each galaxy has a 
matching PSF in the corresponding PSF image 



3. Applying gfit to the GREAT10 data 

3.1. The GREAT10 Galaxy challenge 

We describe in this section the specific pipeline we used in 
the GREAT10 Galaxy Challenge competition that took place 
between December 2010 and September 2011. The aim, content 
and rules of the challenge have been described in the GREAT10 
Handbook ( |Kitching et al.|[20TT| |2012| i. In a nutshell, the main 
goals of the GREAT Challenges are (i) to test existing weak 
lensing measurements methods and (ii) to promote the develop- 
ment of new, more accurate, shear measurement techniques. 

The data consist of 24 datasets of 200 simulated galaxy 
images, each containing 10,000 noisy, PSF-convolved 48 X 48 
pixel galaxy postage stamps, arranged on a 100 x 100 grid (see 
Fig. [T). The GREAT10 edition includes spatially-varying PSF 
and shear fields, contrary to its predecessor, the GREAT08 chal- 
lenge ( |Bridle et al.|2008[ [Bridle et al.|2010) , where these fields 
were set as constant. 

Each of the 24 sets is designed to evaluate the ability of com- 
peting methods to deal with galaxy or PSF fields with different 
properties (e.g. size, signal to noise ratio). We have reproduced 
in Table 2] the main PSF and galaxy characteristic attached to 
each of the GREAT10 set, as specified in the Galaxy Challenge 
results paper (Kitchi ng et al.|2012) , Appendix D. 

The sets were also classified into "Single epoch", "Multi- 
epoch" and "Stable single epoch", depending on whether the in- 
trinsic ellipticities and PSF keep the same or change their spatial 
distribution between images in a set (see Table[T). 

Participating methods were ranked according to the follow- 
ing metrics : 

- A "Raw" quality factor Q according to which the live leader 
board was scored and that measures the difference, averaged 
over all sets, between the reconstructed and true shear power 
spectra. 

- A quality factor Q c {„ obtained after estimation and correction 
of pixel noise. 

- A quality factor Qd„ & train after application of an additional 
training step on top of pixel-denoising. That step consists es- 
sentially in estimating the multiplicative and additive biases 
on high S/N galaxy images (set 7) and applying that calibra- 
tion on the remaining sets. 

- A split of the total bias into an shear-independent additive 
bias and a shear-dependent multiplicative bias between the 
measured and true shear, respectively denoted as c and m 
where c = (c,) and m = (m,), i = (1,2). The c, and m,- are 
respectively the additive and multiplicative biases over the 
two component of shear j\ and ji and are similar to those 
used in STE P I ( |Heymans et al.|2006| > and GREAT08 ( |Bridle| 
|et al.|2 010) for a constant shear field. They are calculated by 



PSF 
Image 001 



PSF 
Image 200 



To each galaxy image corresponds a PSF 
image with 10000 PSF postage stamps 



Fig. 1. The GREAT10 Galaxy challenge dataset structure. There 
are 24 galaxy sets, each containing 200 galaxy images with 
10000 galaxy postage stamps each. To each galaxy postage 
stamp corresponds a matching PSF postage stamp at the same 
spatial position in the galaxy image. 

assuming the error on the shear, y, - j" ue , and the true shear 



yj" ,c obeys a linear relationship of the form 

a = i,2) 



n-yr = m i y>r e +c, 



(3) 



More details are provided in Appendix B of (Kitching et al. 
|2012| i. See e.g. Huterer et al. (2006) for a discussion on the 
origin of additive and multiplicative errors in weak lensing 
studies. 

Additional bias metrics 3i — cr 2 (c) and M - m 2 + 2m, in- 
tended to measure the additive and multiplicative biases cal- 
culated at power spectrum level. Unlike c and m, these met- 
rics account for spatial variability. 



3.2. The GREAT10 gfit implementation 

The GREAT 10 version of gfit only implements the first two steps 
described in Sect. 2.1 that is, PSF correction and galaxy shape 
measurement. 

In GREAT 10, the estimation of the shear field (third step in 
Sect. ~~ 



2.1 



was not mandatory as participants were allowed to 
supply for each image a catalog of estimated galaxy elliptici- 
ties instead of a shear power spectrum (Kitching et al .|201 1) : an 
analysis program was written by the GREAT10 team to calculate 
a shear power spectrum from user-supplied ellipticity catalogs. 
Consequently, like most other competing methods, gfit only pro- 
vided its estimates in the form of a catalog of estimated elliptic- 
ities at requested positions within the images. Future version of 
gfit will allow the extraction of a spatially varying shear. 

The gfit implementation used in GREAT10 consisted of the 
following stages: 

Table 1. Spatial variability of PSF and intrinsic ellipticities: im- 
ages within a set may keep or not the same PSF or intrinsic el- 
lipticity pattern of variation. 



Type of variability within a set 


PSF 


Intrinsic ellipticity 


Type 1 : "Single epoch" 


Variable 


Variable 


Type 2: "Multi-epoch" 


Variable 


Fixed 


Type 3. "Stable single epoch" 


Fixed 


Variable 
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Table 2. Some of the PSF and galaxy properties characterizing 
the GREAT 10 image sets. The second and third columns specify 
whether the PSF or intrinsic ellipticity field were kept constant 
for all images within a set. The parameters in the fourth column 
have been detailed in Kitching et al. ( 2012 1. The default signal to 



Table 3. Main metrics for "gfit den" and "gfit' 



noise (S/N) ratio is 20, while low and high S/N ratios are 10 and 
40 respectively. All sets, except the last four, have galaxies with 
co-centered bulges and disks with a 50/50 bulge-to-disk ratio. 



Set 


PSF 


Intrinsic ellipticity Property of images 


1 


Variable 


Variable 


Fiducial 


2 


Fixed 


Variable 


Fiducial 


3 


Variable 


Fixed 


Fiducial 


4 


Variable 


Variable 


Low S/N 


5 


Fixed 


Variable 


Low S/N 


6 


Variable 


Fixed 


Low S/N 


7 


Variable 


Variable 


High S/N 


8 


Fixed 


Variable 


High S/N 


9 


Variable 


Fixed 


High S/N 


10 


Variable 


Variable 


Smooth S/N 


11 


Fixed 


Variable 


Smooth S/N 


12 


Variable 


Fixed 


Smooth S/N 


13 


Variable 


Variable 


Small galaxy 


14 


Fixed 


Variable 


Small galaxy 


15 


Variable 


Variable 


Large galaxy 


16 


Fixed 


Variable 


Large galaxy 


17 


Variable 


Variable 


Smooth galaxy 


18 


Fixed 


Variable 


Smooth galaxy 


19 


Variable 


Variable 


Kolmogorov PSF 


20 


Fixed 


Variable 


Kolmogorov PSF 


21 


Variable 


Variable 


Uniform bulge/disc ratios 


22 


Fixed 


Variable 


Uniform bulge/disc ratios 


23 


Variable 


Variable 


50/50 bulge/disc offset 


24 


Fixed 


Variable 


50/50 bulge/disc offset 



1. Optional denoising of the galaxy and PSF images with the 



DWT- Wiener method presented in Sect. 2.7 

2. Centroid estimation using SExtractor (Bertin & Arnouts 
|1996| l as described in Sect. |2.6| 

3. Galaxy shape measurement with the gfit program, config- 
ured to use the CCD minimization algorithm described in 
Sect. |2.5| We used 24 x 24 pixel cutouts instead of the full 
48 x 48 pixel original galaxy postage stamps. Similarly, the 
size of PSF postage stamps was reduced to 12 X 12 pixels. 
That decision was made in order to keep the overall compu- 
tation time within acceptable limits and avoid picking-up too 
much noise near the borders of the postage stamps. 

Because of the large number of galaxies, running the pipeline 
on one single processor would not have allowed to meet the 
GREAT10 deadline. Even with a processing time per galaxy of 
0.5 seconds, it would have taken about one month to complete 
process the full GREAT10 dataset. The ability to simultaneously 
run multiple program instances use of parallelism is thus im- 
perative and all programs (denoising, SExtractor wrapper, gfit) 
are written to take advantage of of parallel computers through 
the Message Passing Interface (MPI) ( |Forum| 1995] [1998) . When 
only a few processors are required, the same programs can also 
run on machines with symmetric multiprocessing (SMP) archi- 
tecture. It took about 5 days to process the entire GREAT10 chal- 



Method 


Q 


Qdn 


Qdn & train 


M/2 x 10~ 2 


y/M x 10~ 4 


gfit den 


103.81 


197.88 


229.19 


-2.067 


+0.061 


gfit 


50.11 


122.74 


249.88 


+0.583 


+0.057 



lenge images on a 64-processor machine, which corresponds to 
a processing time between 1 and 2 seconds per galaxy. 

The pipeline is implemented in Python, a programming lan- 
guage known for its power, flexibility and short development 
cycle. The usual standard Python libraries are used, notably: 
NumPy, SciPy, PyFITS and matplotlib. SciPy is the standard 
scientific library for Python and most of its functions consist of 
thin Python wrappers on top of fortran, C and C++ functions. 
SciPy takes advantage of installed optimized libraries such as 
LAPACK (Linear Algebra PACKage) library ( |Anderson et al. 
[1990) . 



4. Analysis of the gfit GREAT10 results 

We summarize and analyze in this section the main Galaxy 
Challenge results as far as gfit is concerned. An overview of the 
GREAT10 results for available participating shear measurement 
methods has already been performed in the GREAT 10 Galaxy 
challenge paper (Kitchi ng et al.|20 12). Our objective here is to 
provide a more detailed analysis of the gfit results. 

We do not, however, analyze the influence of the pixel- 



denoising and training calibration schemes applied in (Kitching 
|et al.|2012) , which we leave for future investigation. 

4.1. Overall results 

The results of the best 12 methods that participated in the 
GREAT 10 Galaxy C hallenge are listed in Table 3 of the 
Kitching et al. ( 2012| l GREAT10 result paper. That list aggre- 
gates results submitted before the official challenge deadline as 
well as submissions made during the so-called "Post challenge", 
a one-week extension to the competition following the deadline. 

Two versions of gfit were submitted during the challenge, 
one named "gfit den cs" that included a denoising step using 



the DWT- Wiener algorithm described in Sect. 2.7 and the other, 
simply named "gfit", that did not. The results obtained by both 
methods are shown in Table [3] 



The gfit version presented in Kitching et al. (20121 is what 
we refer to here as "gfit". In addition, we also present "gfit 
den cs", described in Appendix E5 of Kitching et al. (2012) 
but whose results were not included in the analysis. The "gfit" 
version is identical except that no denoising was applied to 
the data before applying the shape measurement algorithm. To 
simplify, we shorten he name "gfit den cs" to "gfit den" in the 
remainder of this article. 

It can be seen from Table [3] that "gfit den" reaches a raw 
quality factor Q twice as high as that of "gfit". This illustrates 
the gain in accuracy provided by the DWT- Wiener denoising al- 
gorithm. This is further analyzed in Sect. 4.2.2 When the pixel- 



level denoising algorithm of |Kitching et a . ^2012 ) is applied, the 
Qd„ quality factors of both "gfit" and "gfit den" are improved 
by a factor ~ 2, "gfit den" scoring the best Qd n of all methods 
(Qdn = 197.88). The training calibration further increases both 
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"gf it" Q lactors per set (Q) = 50.11 



"gfit den" Q factors per set (Q) = 103.81 



HIM 



t?7CI Fixed PSF variation 

Fixed ellipticity variation 
I 1 Variable PSF/ellipticity 



Fiducial 
Low S/N 
High S/N 
Smooth S/N 
Small Galaxy 
Large Galaxy 
Smooth Galaxy 
Kolrrtogorov PSF 
Uniform b/d 
Offset b/d 




9 10 11 12 13 14 15 16 17 18 19 20 21 22 

Set 




9 10 11 12 13 14 15 16 

Set 



17 18 19 20 21 22 23 24 



Fig. 2. Quality factors per set for "gfit" (left) and "gfit den" (right). The various colors and patterns in the legend indicate the types 
of features simulated in the sets. The acronyms S/N, b/d respectively refer to the signal to noise ratio and galaxy bulge/disc ratio 
or offset. The fiducial S/N was 20 whereas the low S/N and high S/N were respectively set to S/N = 10 and S /N = 40. The 
labels "Fixed PSF" and "Fixed intrinsic ellipticity" correspond to sets where the PSF and intrinsic galaxy ellipticities were spatially 
varying across the field but that variation did not change between images within a set. Further details on the structure of the Galaxy 
challenge dataset and the procedures for calculating the quality factor can be found in in the GREAT10 Galaxy challenge results 
paper Kitching~et al.|(|2012i. 



Qdn&tmin quality factors, especially that of "gfit" (two-fold in- include the "gfit den" power spectrum in Fig. B.2 as well, 
crease). 

We have also included in Table [3] the average additive and 
multiplicative biases 3i and M/2 over all 24 sets. Comparing 

|201 21l result paper, 
additive bias 



with Table 3 of the 
we see that "gfit" reached the 
( yfm = 0.057 x 10~ 4 



Kitching et al.| 

lowest average 



) and the second lowest average multiplica- 
tive bias (At/2 = 0.583 x 10~ 2 ) of all twelve methods (see also 
the plot in Figure. 1, page 6 of that paper). 

We stress that, contrary to what is suggested in Sect. 4.4 of 
|Kitching et al.| ( |2012| l, the low overall bias of "gfit" is intrinsic 
to the method and does not result from the application of a de- 
noising step. The DWT- Wiener algorithm was only used in "gfit 
den", not "gfit" and actually, the average multiplicative bias of 
"gfit den" (M/2 = -2.067 x 10~ 2 ) is higher than that of "gfit". 
Moreover, both methods have similar additive biases. Drawing a 
more refined conclusion about these biases requires an analysis 



at individual set level, which we perform in Sect. 4.3 



4.2. Method accuracy 

In this section we use the quality factor as a measure of accuracy 
and assess the influence of: 

- Galaxy and PSF characteristics included in the images (size, 
signal-to-noise, etc.) 

- Denoising with the DWT- Wiener algorithm 

The quality factors scored for each individual image set are 
plotted on the left-hand side part of Fi g. [2] f or each gfit variant. 



They are also quoted in Tables A. 1 and A. 2 

Due to an editorial mistake, the shear power spectra at- 
tributed to "gfit" in the GREAT10 Galaxy Challenge paper, 
Figure E9, is that from the "fit2 -unfold" method. The correct 
picture was made available at the time of publication and can be 
found here We also provide the correct figure in Fig. |B.1| We 



We focus first on the influence of galaxy and PSF features on 
accuracy. We leave aside the effects of DWT- Wiener denoising 
for now and thus base our analysis on the results of the "gfit" 
variant which is devoid of built-in denoising scheme. The data 
of interest are summarized in Figs.|2]and B.l and in Tables A.l 



4.2.1. Influence of Galaxy and PSF characteristics 

- Influence of signal-to-noise ratio: this is best reflected on 
the results of the "gfit" variant, since it has no built-in noise 
correction scheme. All that was done regarding noise was to 
cut out galaxy and PSF postage stamps to lower dimensions: 
12 x 12 pixels for the PSF and 24 x 24 pixels for galaxies. 
As expected, higher S/N yields higher Q factors: we observe 
a roughly linear progression with a ~ two-fold increase 
from low S/N = 10 to fiducial S /N = 20 and a ~ three-fold 
increase from low S /N — 10 to higher S/N = 40 (see sets 4 
to 12 in Fig. [2}. 



Influence of galaxy size: as seen from the results of sets 
13 to 18 in in Fig.|2]i, "gfit" seems quite sensitive to galaxy 
size, the worst Q factors being obtained on smaller galaxies 
and the best on the larger ones, with a factor ~ 2 difference. 
The use of postage stamp cut-outs of identical dimensions, 
regardless of the actual FWHM of the galaxy they contain 
could be responsible for this effect, as cut-outs with smaller 
galaxies are likely to be more noise-dominated than those 
with larger objects. It may also be that the minimizer is 
less accurate on smaller objects. "Smooth" galaxies have 
sizes varying according to a Rayleigh distribution ( Kitching 
et al.|2012) , so it is not surprising that the Q factors of the 
corresponding sets 10, 11, 12 take values in between those 
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of small and large galaxies. 



4.3. Bias analysis 



Influence of bulge/disc distribution and offset: having 
varying b/d ratios (sets 21, 22) seems to decrease accuracy 
down to the level of small-size galaxies (sets 13, 14). On the 
other hand, introducing a non-zero 50/50 bulge-to-disc (b/d) 
offset (sets 21, 22) tends to yield slightly higher accuracy 
compared to "smooth" galaxies (sets 17, 18). So it seems 
the single-component galaxy model (see Sect. 2.3 i of "gfit" 
is more sensitive to b/d ratio than to b/d separation. The use 
of a more sophisticated galaxy model would certainly yield 
a small gain for some types of galaxies but at the price of 
an additional computational time for model fitting. On real 
data, the importance of the b/d distribution and offset will 
probably depend on the available galaxy sample. 

Influence of turbulence: The raw Q factor plot in Fig. [2] 
does not show a strong impact from the inclusion of a 
Kolmogorov power spectrum in the PSF ellipticities. We 
note however a higher score on the fiducial turbulent set 
compared to the non-turbulent one for no obvious reason. 
This seems counter intuitive as PSF turbulence usually 
degrades accuracy and this phenomenon may actually not 
be related to turbulence. 

Influence of spatial variability between images: as de- 
scribed in Kitching et al.| ( |2012| l, the Galaxy challenge data 
are divided into so-called "Single epoch", "Multi-epoch" and 
"Stable single epoch" depending on whether the intrinsic el- 
lipticities and PSF keep the same or change their spatial dis- 
tribution between images in a set (see Table[T]). It seems that 
having fixed instead of variable intrinsic ellipticities ("Multi- 
epoch" sets) slightly decreases accuracy. Apart from this, no 
clear trend really stands out from the results and it is not 
clear whether the difference in accuracy is due to the type 
of spatial variability used or to the specific sample of im- 
ages chosen for a set. This topic nevertheless deserves to be 
investigated further in a separate work. 



4.2.2. Effect of denoising on accuracy 



We discuss in this section the effect on accuracy of the ap- 
plication of the DWT- Wiener denoising scheme described in 
Sect. |2.7| The influence of denoising can be clearly observed by 
comparing the plot of "gfit den" (left) with that of "gfit" in F ig. [2 
(right). The corresponding scores are also listed in Tables A.l 
andlA~2l 



We find an average two-fold increase in accuracy, the effect 
being stronger on high S/N images and larger galaxies. Sets with 
small galaxies are only slightly improved, however. We also no- 
tice that the plots of "gfit" and "gfit den" show identical quality 
factors for low S/N sets 4 to 6. Further investigation showed that 
denoising was, by mistake, not applied on those sets. This would 
probably have improved the overall Q factor of "gfit den". 

These results strongly suggest that the DWT- Wiener algo- 
rithm really improves the overall accuracy on galaxy shape mea- 
surement. We also note that denoising does not alter the quality 
factor hierarchy between sets: the sets with best scores in the 
"gfit" plots remain the same in the "gfit den" plot. 



We investigate in this section how multiplicative and additive 
biases are affected by galaxy properties and the use of denoising. 

4.3.1. Influence of Galaxy and PSF characteristics 

To complement the results of Table [3] relative to bias, we have 
plotted in Figs. [3] and [4] the multiplicative and additive biases of 
each set. As noted in Sec. |4.1| "gfit den" reached the lowest av- 
erage additive bias and the second lowest average multiplicative 
bias of all 12 twelve competing methods. 

Focusing on the multiplicative bias and leaving aside the ef- 
fect of denoising for now, we can make a few observations from 
the left-hand side plot of Figs. [3] 

- Highest and lowest multiplicative bias: all multiplicative 
bias values, except one (set 5), lie below M/2 = 3.5 X 10" 2 
in absolute value. The largest biases are found on fiducial set 
1 (fixed PSF & variable ellipticities), set 5 (low S/N, fixed 
PSF) and set 18 (smooth galaxies, fixed PSF). The smallest 
biases are obtained on set 4 (low S/N, fixed PSF & variable 
ellipticities), set 17 (smooth galaxies, fixed PSF) and sets 
22 to 24 (galaxies with varying b/d ratios and fixed PSF, 
galaxies with non-zero b/d offset). 

- Influence of signal-to-noise ratio: the results from sets 1 
and 7 suggest, as expected, that multiplicative bias decreases 
with S/N. The low S/N set 4, however, shows a very small 
bias. It may be that the true bias is large and positive but that 
it was offset by e.g. a negative bias due to the contribution 
of large galaxies in that particular set. The relatively large 
bias of set 5 compared to sets 4 and 6 may also be an 
consequence of the PSF having a fixed variation pattern in 
that set. 

- Influence of galaxy size: small galaxies (sets 13, 14) 
and large galaxies (sets 15, 16) have comparable biases 
in absolute value, smaller galaxies having a positive bias 
and larger ones a negative bias. A lower bias is reached 
on "Smooth" galaxy images (set 17), likely because the 
negative and positive biases of small and large galaxies 
respectively compensate each others. The bias in set 18 may 
also have been artificially amplified by the fixed variation 
pattern of the PSF in that set. 

- Influence of bulge/disc distribution and offset: letting 
the b/d ratios of galaxy vary within a set tends to yield a 
multiplicative bias comparable to that of small-size galaxies, 
probably causing the decrease in accuracy mentioned in 
Sect. |4.2| In contrast, the biases associated with a non-zero 
50/50 b/d offset are among the lowest. As noted earlier, the 
underlying single-component Sersic-based galaxy model of 
"gfit" (see Sect. 2.3 i seems handle quite well profiles with 



an off-centered bulge and disk. 

Influence of turbulence: the introduction of PSF with 
turbulent ellipticities (set 19) does not induce a particularly 
significant bias compared to the average. 

Influence of spatial variability between images: we note 
that having "Fixed PSF" and "Fixed intrinsic ellipticity" sig- 
nificantly alters the multiplicative bias, especially on low 
S/N and smooth galaxy sets. Because all images have the 
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Fig. 3. Multiplicative bias M/2 per set. The legend patterns and colors are identical to those of Fig. [2] 



xio- 5 "gfit" additive bias per set (v^) = 5.73 x 10~ s x io-« "gfit den" additive bias per set (\/ r A) = 6.15 x 10" 15 




Fig. 4. Additive bias V3^ per set. The legend patterns and colors are identical to those of Fig. U. 



same spatial variation within a set, its is likely that the bias 
of one image is just amplified as many times as there are 
images (i.e. 200 times). The effect is particularly strong on 
image with low S/N (set 5) and smooth galaxies (set 18). 

As far as additive bias is concerned, we note the following trend: 

- Highest and lowest additive bias: "gfit" obtains an additive 
bias < 10~ 5 on all types of sets, except on the low 
S/N ones. All values are positive. The lowest bias is reached 
on small galaxies (sets 13, 14) and the largest on low S/N 
images (sets 4 to 6). 

- Influence of signal-to-noise ratio: as for multiplicative 
bias, we find a trend toward higher biases for lower S/N 
(sets 1 to 9). 

- Influence of turbulence: the introduction of a Kolmogorov 
power spectrum in PSF ellipticity induces a slightly greater 
additive bias than average (see set 19). 



- Influence of galaxy size: as for multiplicative bias, we note 
a tendency of smaller galaxies to have a higher additive bias, 
as reflected by the values for sets 15 to 18. 

- Influence of bulge/disc distribution and offset: the additive 
bias appears larger on galaxies with varying b/d ratios (sets 
21, 22) and non-zero b/d offset (sets 23, 24). 

4.3.2. Effect of denoising on bias 

Comparing the plots for "gfit" and "gfit den" in Figs.|3]and|4] we 
find that the additive bias does not change significantly, keeping 
about the same bias values per set. In contrast, the structure of 
the multiplicative plots is significantly altered. 

As seen on the "gfit den" plot, denoising tend to intro- 
duce some amount of negative multiplicative bias on all sets. 
Although the amount of bias on fiducial sets remains roughly 
the same in absolute value, the DWT- Wiener algorithm clearly 
impacts the multiplicative bias relative to galaxy size, b/d ratio 
and turbulence. 
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As regards galaxy size, even though the multiplicative bias 
on small galaxies is almost unchanged, that on large galaxies 
increases about six-fold. Because of that, the bias on "Smooth" 
galaxies, which also includes large galaxies, also increases. 

DWT- Wiener denoising seems to also improves the resolu- 
tion of the bulge and disk components, so that "gfit" has more 
difficulty fitting its underlying single-component Sersic model to 
images with "uniform" and "offset" b/d. The effect is stronger on 
larger galaxies, causing a > tenfold multiplicative bias increase. 

Lastly, the introduction of PSF turbulence result in a ~ five- 
fold multiplicative bias degradation. 

Despite the degradation of multiplicative bias on some 
sets, the accuracy of shape measurements increases two-fold as 
shown by the corresponding gain in Q factor. We also note that 
denoising improves the results on sets that already have a high 
S/N. All in all, the use of denoising is thus clearly beneficial. 

5. Conclusions 

We have described in this paper the gfit shape measurement 
method, a model-fitting based on a simple Sersic galaxy model 
(Sect. [2]). The method uses a custom-developed minimizer based 
on a "coordinate descent" algorithm that finds a local minimum 
with the lowest x 1 of the residuals between true and modeled 
galaxy. 

We have also performed an analysis of our results in the 
GREAT 10 Galaxy Challenge (Sect. [3). We participated in the 
competition with two gfit variants: "gfit den", which applied 
a denoising step before performing model-fitting, and "gfit", 
which did not use denoising. The noise removal technique 
employed is DWT-Wiener, a wavelet-based, shape-preserving 
algorithm particularly suitable for shape measurement, (see 
Sect.|T7). 

We highlight below the main conclusions of your analysis. 

- Method accuracy: accuracy improves significantly as 
S/N gets higher and galaxy size larger. The underlying 
simple Sersic-based galaxy model of gfit has more difficulty 
handling galaxies with a non-zero 50/50 offset between 
bulge and disc. The inclusion of Kolomogorov turbulence 
in ellipticities is not seen to yield a significant change in 
accuracy. 

- Additive and multiplicative bias: the non-denoised "gfit" 
variant reached the lowest average additive bias and second 
lowest average multiplicative bias of all twelve competing 
methods. Both additive and multiplicative bias tend to be 
larger on galaxies with high S/N, smaller size and galaxies 
with b/d ratio differing from 50/50. 

- Impact of denoising: the application of the DWT-Wiener 
noise removal algorithm yields a two-fold improvement in 
accuracy (Q factor) despite significantly degrading the mul- 
tiplicative bias on galaxies with a high S/N, small size and 
significant bulge/disk ratio and separation. 

It is interesting to see that, despite the simplicity of the 
galaxy model used, its results in the Galaxy Challenge estab- 
lished gfit as one of the four top-performing methods, both in 
terms of accuracy and bias. Given that the results were obtained 
on simulated data, this raises the question of how important 
having a realistic galaxy model really matters when measuring 



galaxy shapes from real data. Providing more clues on this 
question is one of the objectives of the forthcoming GREAT3 
challenge. 
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Table A.l. "gfit": results per set. Sets with S/N 10, 40 are re- 
spectively highlighted in orange and blue. Fiducial sets with S/N 
20 are represented in green and all remaining sets also have a 
S/N of 20. 



SCt 
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Qdn 


*Zdn & train 


A A /O w 1 C\— 2 

JWI 1 x to 


yJx x 10 
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+0.054 
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-0.432 


+0.103 
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+0.150 
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31.32 
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-2.136 


+0.118 
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+2.270 


+0.041 
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120.18 
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+0.176 


+0.028 
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185.17 


184.27 


+ 1.780 


+0.011 
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+0.057 


1 1 
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17 


52.88 


125.30 


246.64 


+0.334 
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Table A.2. "gfit den": results per set. Sets with S/N 10, 40 are 
respectively highlighted in orange and blue. Fiducial sets with 
S/N 20 are represented in green and all remaining sets also have 
a S/N of 20. 
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59.06 
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+0.671 
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198.34 
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158.18 
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All 


103.81 


197.88 


229.19 
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Appendix A: Accuracy and bias per set 



Tables [ATT] and A. 2 respectively quote the actual quality factor 
and bias values reached by the non-denoised and denoised vari- 
ants of the gfit shape measurement method. 



Appendix B: Sear power spectra 



Figs |B.l| and B.2 respectively show the shear power spectra of 
the non-denoised and denoised variants of gfit submitted in the 
GREAT 1 galaxy Challenge . 
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Fig. B.l. The shear power spectra of the "gfit" variant over the 24 sets of the Galaxy challenge. The red lines denote the estimated 
shear power spectra while the green lines represents the true shear power. The blue lines indicate the spectra after application of 



pixel-level denoising (not DWT- Wiener denoising). The code run to plot these power spectra is identical to that used in the (Kitching 
et al.|2012 1 paper. Note that the power spectra attributed to "gfit" in the Figure E9 of ( |Kitching et al.|2012] > are actually those of the 
"fit2-unfold" method. 
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Fig. B.2. The shear power spectra of the "gfit den" variant over the 24 sets of the Galaxy challenge. This variant has a preliminary 
denoising step using the DWT- Wiener algorithm. The red lines denote the estimated shear power spectra while the green lines 
represents the true shear power. The blue lines indicate the spectra after application of pixel-level denoising (not DWT- Wiener 
denoising). The code run to plot these power spectra is identical to that used in the (Kitching et al.|2012 1 paper. 
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